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1.  INTRODUCTION 
A.  THE  PROBLEM 

An  efficient  procedure  is  developed  in  this  report  for  analyzing 
scattering  of  an  electromagnetic-plane  wave  from  a  homogeneously  loaded 
rectangular  cavity  recessed  in  a  perfectly  conducting  ground  plane.  The 
analysis  is  performed  for  the  two-dimensional  geometry  shown  in  Figure 
1.1.  The  interest  in  this  work  stems  in  large  part  from  the  significant 
radar  cross  section  (RCS)  contribution  of  such  geometries.  Of  special 
interest  are  loaded  and  unloaded  notch  geometries,  in  which  the  width 
and  depth  of  the  rectangular  cavity  are  small  to  moderately  large  in 
terms  of  the  wavelength,  and  also  wide  shallow  cavities  in  which  the 


Figure  1.1.  Electromagnetic  fields  incident  on  a  loaded  cavity  recessed 
in  a  ground  plane. 
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width  of  the  cavity  is  large  in  tenns  of  the  wavelength  while  the  depth 
is  a  fraction  of  a  wavelength.  The  latter  geometry  is  of  particular 
interest  because  of  its  its  application  to  surface  wave  antennas,  and  to 
arrays  of  dielectric  covered  slots. 

The  notch  geometry  is  analyzed  by  a  conventional  numerical  moment 
method  technique  cast  into  a  particular  efficient  form.  The  wide 
shallow  cavity,  because  of  its  large  electrical  size,  cannot  be 
efficiently  analyzed  by  the  same  technique.  Thus,  a  hybrid  version  of 
the  moment  method  is  developed  for  this  geometry  which  results  in  a 
considerable  savings  of  CPU  time  compared  to  the  conventional  moment 
method. 


B.  PREVIOUS  RESEARCH 

We  consider  first  some  previous  work  done  on  problems  related  to 
the  scattering  of  a  plane  wave  from  2-D  dielectric  loaded  notch 
geometries.  In  each  of  these  solutions  the  moment  method  [1]  was 
employed  to  solve  an  integral  equation  involving  the  equivalent  magnetic 
current  in  the  aperture.  Because  of  the  relatively  small  electrical 
size  of  the  notches  conventional  moment  method  techniques  were  used. 
Richmond  [2]  solved  the  problem  of  a  dielectric  loaded  notch  recessed  in 
a  perfectly  conducting  circular  cylinder  by  expanding  the  fields  in  the 
dielectric  in  a  set  of  cylindrical  modes  and  then  applying  Galerkin's 
method  to  determine  the  aperture  fields.  Kautz,  Pathak  and  Peters  [3], 
as  well  as  Wang  [4]  considered  the  geometry  of  a  dielectric  loaded  gap 
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in  a  thick  perfectly  conducting  screen.  They  determined  the  equivalent 
magnetic  current  (in  both  the  upper  and  lower  apertures)  by  employing 
the  moment  method  in  conjunction  with  the  multiple  scattering  method. 

By  placing  an  impedance  wall  in  the  lower  aperture  and  assigning  to  it 
an  impedance  of  zero,  they  are  able  to  treat  a  loaded  notch  geometry* 
However,  this  method  of  solution  for  the  notch  problem  is  not  very 
efficient,  primarily  because  of  the  added  complexity  resulting  from  the 
formulation  which  allows  one  to  treat  the  more  general  case  of  loaded 
gaps  in  thick  screens. 

We  now  consider  some  relevant  previous  research  done  on  problems 
related  to  loaded  wide  shallow  cavities.  Pathak  and  Kouyoumjian  [51 
solved  the  canonical  problem  of  surface  wave  diffraction  by  a  truncated 
dielectric  slab  recessed  in  a  perfectly  conducting  gound  plane.  They 
restricted  their  analysis  to  a  lossless  dielectric  slab  able  to  support 
only  the  TMq  surface  wave.  The  reflected  and  diffracted  fields  are 
expressed  in  terms  of  the  geometrical  theory  of  diffraction  (RTD)  [61. 
The  surface  wave  reflection  and  diffraction  coefficients  are  obtained 
from  a  formally  exact  solution  employing  a  combination  of  the 
generalized  scattering  matrix  technique  (GSMT)  [71  and  the  Wiener-Hopf 
procedure  [81.  This  solution,  although  formally  exact,  requires 
considerable  computation.  A  more  efficient  solution,  which  can  also 
handle  a  discretely  inhomogeneous  portion  of  dielectric  near  the 
truncation,  is  given  by  Chuang  [81.  This  solution  employs  a  moment 
method  solution  to  an  integral  equation  involving  the  equivalent 
magnetic  current  at  the  air-dielectric  interface.  Because  of  the 
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semi -infinite  geometry,  it  is  necessary  to  use  a  hybrid  moment  method 
[lOl,  [111  solution  in  which  the  basis  functions  are  pulses  near  the 
truncation  and  in  the  inhomogeneous  portion  of  the  dielectric;  whereas, 
away  from  the  truncation  and  inhomogeneous  regions,  the  unknown  surface 
field  is  expressed  in  terms  of  the  incident  and  reflected  surface  waves 
together  with  the  diffracted  wave  (which  vanishes  to  first  order  on  the 
dielectric  interface).  The  forms  of  the  reflected  and  diffracted  waves 
are  known  and  their  amplitudes  constitute  unknowns  to  be  determined.  We 
can  use  this  solution  in  conjunction  with  reciprocity  to  determine  the 
lauching  coefficient  of  a  surface  wave  due  to  a  plane  wave  incident  on 
the  canonical  geometry.  However,  in  this  situation,  reciprocity  will 
not  yield  any  information  on  the  surface  transition  region  which  extends 
a  couple  of  wavelengths  (and  even  farther  for  a  plane  wave  incident  at 
grazing  angles)  from  the  truncation.  This,  coupled  with  the  fact  that 
the  semi -infinite  portion  of  dielectric  must  be  lossless,  limits  the 
usefulness  of  that  canonical  problem  in  obtaining  a  solution  to  the 
problem  shown  in  Figure  1.1.  In  addition,  the  solutions  of  Pathak,  et 
al.  and  Chuang  are  only  for  an  incident  TMq  surface  wave;  the  other 
polarization  is  not  treated,  simply  because  the  TEj  surface  wave  is 
cutoff  for  the  configuration  considered  by  Pathak,  et  al .  in  [51. 

Another  solution  relevant  to  the  loaded  shallow  cavity  problem  is 
given  by  Richmond  [12]  for  scattering  of  a  plane  wave  from  thin 
dielectric  strips  of  infinite  length;  the  electric  field  is  parallel  to 
the  edges  of  the  strip.  The  field  in  the  dielectric  is  expanded  in  a 
sum  of  forced  wave  and  two  surface  waves,  with  the  amplitudes  of  these 
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waves  being  determined  by  Galerkin's  method.  This  solution  Is  relevant 
to  the  loaded  shallow  cavity  problem  because  of  the  physical  nature  of 
the  basis  functions.  However,  the  presence  of  the  perfectly  conducting 
surface  In  the  case  of  the  shallow  cavity  yields  diffraction  effects 
that  are  significantly  different  from  the  diffraction  characteristics  of 
the  dielectric  strip. 

Finally,  we  mention  a  hybrid  moment  method  solution  to  a  problem 
not  related  to  that  being  considered  in  this  report,  but  which  well 
illustrates  the  computational  capabilities  of  the  technique.  Srikanth, 
Pathak,  and  Chuang  [131  employed  a  hybrid  uniform  theory  of  diffraction 
(UTD)  [14],  [151,  -moment  method  (MM)  solution  to  analyze  the  plane  wave 
scattering  from  an  infinite,  perfectly  conducting  semi-circular 
cylinder.  A  mixed  basis  expansion  involving  pulses  near  the  edges  and 
UTD  functions  elsewhere  is  used  for  the  unknown  current  on  the  cylinder. 
The  basis  function  amplitudes  are  then  determined  by  point-matching. 
Numerical  results  show  excellent  agreement  with  conventional  pulse  basis 
expansions,  and  are  obtained  in  a  fraction  of  the  CPU  time  of  the  latter 
method.  We  point  out,  however,  that  the  hybrid  UTD-MM  is  most 
applicable  to  a  particular  geometry  when  the  diffraction  characteristics 
(i.e.,  UTD  diffraction  effects)  of  the  canonical  features  of  that 
geometry  are  well  understood. 

In  this  thesis  an  efficient  conventional  moment  method  solution  is 
presented  for  the  scattering  of  a  plane  wave  from  a  notch  geometry  and  a 
very  efficient  hybrid  moment  method  solution  is  presented  for  the 
special  case  of  loaded  wide  shallow  cavity.  For  the  notch  problem  the 
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equivalent  magnetic  current  in  the  aperture  is  expanded  in  a  set  of 
parallel -plate  waveguide  modes  in  which  each  mode  is  an  entire  domain 
basis  function  across  the  aperture.  The  unknown  field  is  only  over  the 
extent  of  the  aperture  because  the  Green's  function  chosen  in  the 
formulation  of  the  integral  equation  for  the  unknown  in  this  geometry  is 
comprised  of  two  terms,  one  for  a  half-space  with  a  perfectly  conducting 
ground  plane,  and  the  other  for  a  loaded,  fully-enclosed,  perfectly 
conducting  cavity.  Galerkin's  method  is  used  to  determine  the  mode 
amplitudes.  Only  3-4  modes  per  free-space  wavelength  are  needed  for  a 
converged  solution,  as  compared  to  the  10-20  pulses  or  piecewise 
sinusoids  per  wavelength  that  are  usually  necessary  if  subsectional 
basis  functions  were  used.  By  using  an  appropriate  form  of  the  cavity 
Green's  function,  the  cavity  contribution  to  an  individual  impedance 
matrix  element  is  arranged  into  a  closed-form  term  plus  a  highly 
convergent  summation.  In  fact,  when  using  the  waveguide  mode  basis 
functions,  the  closed-form  contribution  alone  suffices.  This,  in 
addition  to  an  evaluation  of  the  half-space  contribution  employed  by 
Richmond  [12],  yields  a  very  efficient  matrix  fill.  Excellent  agreement 
has  been  found  between  the  modal  basis  expansion  presented  here  and  a 
standard  pulse  basis  solution.  In  addition,  the  modal  basis  solution 
was  observed  to  be  several  times  faster  than  the  pulse  basis  solution. 
Practically  speaking,  the  modal  basis  solution  represents  a 
cost-effective  solution  (on  a  VAX  11/780)  to  scattering  from  recessed 
rectangular  cavities  of  width  <  25  wavelengths  and  of  arbitrary  depth. 
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Because  of  the  special  interest  in  dielectric  covered  antenna 
cavities,  a  separate  solution  was  performed  for  the  wide  and  shallow 
loaded  cavity  geometry.  By  'wide  and  shallow'  it  is  meant  that  the 
cavity  width  is  greater  than  several  wavelengths  and  the  cavity  depth 
and  loading  are  such  that  at  most  only  one  surface  wave  mode  exists. 

The  same  integral  equation  and  Green's  functions  are  used  as  for  the 
previous  discussion,  but  a  different  set  of  basis  functions  is  employed. 
In  this  hybrid  moment  method  solution  the  equivalent  magnetic  current  in 
the  cavity  aperture  is  expanded  as  a  forced  wave  and  two  surface  waves 
(as  in  [121)  plus  several  additional  basis  functions  that  describe 
transition  effects  which  are  most  significant  near  the  cavity  edges. 
Again,  Galerkin's  method  is  employed,  but  now  the  number  of  unknowns 
remains  almost  constant  (<  10)  regardless  of  the  cavity  width.  Because 
the  transition  effects  are  described  via  a  choice  of  exponential 
functions,  the  computation  of  the  impedance  matrix  is  not  significantly 
different  from  that  of  the  modal  basis  solution.  The  result  is  an 
extremely  efficient  solution  for  the  wide  and  shallow  loaded  cavity. 

Good  backscatter  results  have  been  obtained  for  the  magnetic  field 
parallel  to  the  cavity  edges  for  which  case  the  TM^  surface  wave  always 
exists,  and  when  the  cavity  depth  is  nearly  equal  to  or  less  than  0.1 
wavelengths.  Tor  the  other  polarization,  i.e.,  when  the  electric  field 
parallel  to  the  cavity  edges,  the  lack  of  a  surface  wave  mode  in  the 
dielectric  for  shallow  cavities  causes  the  hybrid  moment  method  solution 
to  be  less  accurate  at  grazing  angles  than  the  corresponding  solution 
for  the  TM  polarization.  For  both  polarizations  the  hybrid  basis 
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solution  is  far  more  computationally  efficient  for  wide  apertures  than 
is  the  modal  basis  solution. 

C.  FORMAT 

In  Chapter  II,  the  integral  equations  involving  the  equivalent 
magnetic  current  in  the  aperture  are  derived.  In  Chapter  III,  the 
moment  method  solution  to  the  integral  equations  is  described.  In 
Chapter  IV,  numerical  results  for  specific  examples  are  presented  and 
discussed.  In  Chapter  V,  the  conclusions  are  stated. 

A  time  convention  of  exp(j<Dt)  is  assumed  and  suppressed;  only  here 
does  t  refer  to  time. 
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II.  INTEGRAL  EQUATIONS  FOR  THE  APERTURE  FIELDS 

A.  LOADED  CAVITY  RECESSED  IN  A  PERFECTLY  CONDUCTING  GROUND  PLANE 

Consider  electromagnetic  fields  incident  on  the  2-D  geometry  shown 
in  Figure  2.1.  The  source  of  the  incident  fields  is  assumed  to  be  in 
region  I  as  shown  in  the  figure.  The  material  in  region  II  is  assumed 
to  be  linear,  homcgeneous,  isotropic  and  time-invariant  and  may  be 
either  lossy  or  lossless.  The  problem  considered  here  is  the 
determination  of  the  fields  scattered  in  region  I  due  to  a  known  source 
which  is  also  in  region  I. 


REGION  I 


Figure  2.1.  An  arbtrary  loaded  cavity  recessed  in  a  perfectly 
conducting  ground  plane. 
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Since  the  geometry  is  2-0,  the  general  vector  problem  can  be  easily 
reduced  to  2  scalar  problems.  The  incident  fields  may  be  decomposed 
into  2  polarizations:  (1)  the  electric  field  vector  parallel  to  the 
conducting  edge,  denoted  as  the  'soft  polarization'  (subscript  's') 
case,  and  (2)  the  magnetic  field  vector  parallel  to  the  conducting  edge, 
denoted  as  the  'hard  polarization'  (subscript  'h')  case.  Thus,  we 
define  the  total  scalar  fields  to  be 


U5  (x,y) 
h 


(x,y) 

(x.y) 


(2.1) 


Integral  equations,  in  which  the  unknown  happens  to  be  the  electric 
fields  only  within  the  aperture  at  y=0  and  |x|  <  a,  may  be  derived  in  a 
straight-forward  manner  by  the  application  of  Green's  second  identity 
[161.  The  details  of  this  derivation  are  deferred  to  Appendix  A;  here 
we  present  only  the  resulting  equations  and  their  interpretation  in 
terms  of  the  equivalence  principle.  These  integral  equations  are  later 
solved  by  a  moment  method  procedure  [11  using  a  special  set  of  only  a 
few  select  basis  functions  as  discussed  in  Chapter  III. 

For  the  total  fields  in  region  I,  we  employ  the  equivalence 
principle  [16]  to  replace  the  aperture  by  an  equivalent  magnetic  surface 
current  source  distribution  just  above  a  perfectly  conducting  surface 
which  short-circuits  the  aperture.  Then,  convolving  this  equivalent 
source  distribution  with  respect  to  a  perfectly  conducting  half-space 
Green's  function,  j/2  (k|p-p'|),  we  find  the  fields  in  region  I 

(y>0)  to  be 
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(x.y)  =  (x,y)  +  (x.y) 


+  /  dx' 
-a 


^  [j/2  H^^^(k|5-p'|)]  .  u^Cx'.y') 
(?)  - 

J/2  (k|p-p‘l)  •  IT-  (x'.y') 


(2.2) 


y'=n 


1  r 

where  U5’  are  the  incident  and  reflected  fields  of  the  equivalent 
h 

(short-circuited)  problem,  and  where 


Ip-p'I  =  •Tx-x')^  +  (y-y')2 


(2.3) 


For  the  fields  in  region  II  (cavity  region)  we  employ  the 

equivalence  principle  to  replace  the  aperture  by  an  equivalent  magnetic 

surface  current  source  distribution  just  below  the  short-circuited 

aperture.  Then,  integrating  the  equivalent  sources  over  the  Green's 

function,  65,  for  the  short-circuited  cavity,  we  find  the  fields  in 
h 

region  II  (y<0)  to  he 


a 

u  (x.y)  =  +  /  dx 
h  -a 


”3Gs 

17-  (x,y,;x' ,y' )  •  u^  (x’ ,y' ) 
9Uh 

_G^  (x,y;x',y')  •  (x'.y')^ 


.  (2.4) 

y'=n- 


The  cavity  Green's  functions,  G5,  are  discussed  in  the  next  section. 

h 

The  equivalent  source  distributions  of  Equations  (2.2)  and  (2.4) 
may  be  readily  interpreted  as  equivalent  magnetic  currents  by  noting 


that  ^5(x'.y')  <_Q  is  the  z-directed  aperture  electric  field  for  the 

3Uh  I 

TM^  or  soft  case  and  that  3y^(x' .y ' ) i.g  is  proportional  to  the 
x-directed  aperture  electric  field  for  the  TE^  or  hard  case. 


U 


respectively.  The  equivalent  magnetic  surface  currents  are  then 
obtained  from 

=  E  X  n  (2.5) 

A  A  A 

where  n  =  y  for  region  I  and  n  =  -y  for  region  II.  The  equivalent 
currents  for  the  2  regions,  denoted  by  and  are  qualitatively 
shown  in  Figures  2.2a  and  2.2b  for  the  hard  and  soft  polarizations, 
respectively.  Note  that  for  the  hard  polarization  the  current  is 
z-directed,  while  for  the  soft  polarization  it  is  x-directed. 


(a)  Equivalent  magnetic  currents  for  hard  polarization. 


(b)  Equivalent  magnetic  currents  for  soft  polarization. 


Figure  2.2.  Equivalent  magnetic  current  sources  in  the  aperture  region 
of  a  cavity  recessed  in  a  ground  plane. 
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To  arrive  at  integral  equations  in  which  the  unknowns  are  the 
aperture  electric  fields,  we  equate  the  fields  as  given  by  Equations 
{2,2)  and  (2.4)  in  the  aperture.  Continuity  of  tangential  E  and  H 
fields  across  the  aperture  implies  the  following  relations; 

u^(x,y=0‘^)  =  u 


3uh 

■5y(x.y) 


1  3Uh 

T;:  37'(x.y) 


(2.6b) 


9Uc 


37"(x,y) 


1 


y=0  ■  Up  3y 


(x.y) 


y=0- 


(2.6c) 


Since  we  ultimately  will  calcuate  the  fields  in  region  I  (y>0),  the 
resulting  integral  equations  are  expressed  in  terms  of  the  aperturt 
fields  at  y^O*.  For  the  hard  polarization  we  obtain  the  integral 
equation 


-a 


/  dx'  [j/2  H^^^(k|x-x' I )  -  G^(x,x’)]  dy 

=  -2  uJ(x.O) 


3uh(x' ,y’ ) 


y*=0 


(2.7) 


Note  that  in  deriving  Equation  (2.7)  we  have  used  Equations  (2.6a) 
and  (2.6b)  to  enforce  the  continuity  of  tangential  H  and  tangential  E, 
across  the  aperture,  respectively. 

To  derive  an  integral  equation  for  the  soft  polarization  we  take  3y 
of  Equations  (2.2)  and  (2.4)  and  equate  the  results  in  the  aperture. 

The  continuity  relations  (2.6a)  and  (2.6c)  are  applied  to  ensure 
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continuity  of  tangential  E  and  tangential  H,  respectively.  The 


/ 

"d 


3  integral  equation  is 

1  32  G5(x,y;x',y') 

3y3y' 

y*=0+  Pr  3y3y' 

y’=o- 

— 

y=0++ 

1 

y=0-- 

^  I  i  r 

Uj(x',0)  dx'  =  -3y  I  u^(x,y)  +  Uj(x,y) 


y=0 


(2.8) 


The  derivative  ^  was  necessary  for  the  enforcement  of  continuity  of 
tangential  H. 

Once  Equations  (2.7)  and  (2.8)  have  been  solved  for  their 
respective  equivalent  source  distributions,  the  scattered  fields  in 
region  I  may  be  determined  from  Equation  (2.2)  with  the  short-circuit 
terms  left  out.  Since  we  are  primarily  interested  in  the  far-zone 
fields,  the  asymptotic  approximation  to  the  HankeT  function  [171 


.-jkp 

TT" 


is  used  along  with  the  usual  approximate  form  of  |p-p'|  [161 


p-p 


p  -  p*p'  in  phase 
p  in  amplitude 


The  scattered  fields  are  then 

(p,(i>)  =  “ 

h 


e  ''  '  “jksin(|>~ 
_  1  _ 


-jkp  ^ 

f _  Jkx'cos^  '  " 


-a 


I  ^u^,(x') 

IJy"  . 


dx' 

(2.9) 
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B.  LOADED  RECTANGULAR  CAVITY  RECESSED  IN  A  PERFECTLY 
CONDUCTING  GROUND  PLANE 

In  this  section  we  consider  the  specific  problem  to  be  solved; 
namely,  the  scattering  of  a  plane  wave  from  a  loaded  rectangular  cavity 
recessed  in  a  perfectly-conducting  ground  plane,  as  shown  in  Figure  2.3. 
For  plane  wave  incidence,  the  incident  fields  may  be  written  in  the 

form 

u!(x.y)  =  1“  .  (J.lOa) 

h  ^0 

Then  for  the  equivalent  short-circuited  problem,  the  incident  fields  are 
still  given  by  Equation  (2.10a),  whereas  the  fields  reflected  from  the 
short-circuit  are  given  by 


Figure  2.3.  Plane  wave  incidence  on  a  loaded  rectangular  cavity 
recessed  in  a  ground  plane. 
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A  rigorous  derivation  of  the  Green's  function  for  the 


short-circuited  rectangular  cavity  is  given  in  Appendix  B.  Here  we 

present  the  governing  relations  and  the  explicit  form  of  for  the 

h 

rectangular  cavity.  The  Green’s  function  must  satisfy  the  differential 
equation 


+  ^2  e^u^]  Gs(x,y;  x',y’)  =  -6(x-x')  6(y-y')  (2.11a) 

n 

for  all  (x.y)  and  (x',y')  within  the  cavity,  and  the  boundary  conditions 

3G^  =  0  on  C  (2. 11b) 

W 


where  C  is  the  contour  which  defines  the  enclosed  cavity  and  n  is  the 
inward  normal  on  C.  For  the  geometry  of  Figure  2.3  the  cavity  is 
defined  by  y=-b,  y=0  and  x*-a,  x=a. 

We  solve  Equation  (2.11a)  by  separation  of  variables  and  an 
application  of  the  boundary  condition  (2.11b).  Note  that  the  boundary 
conditions  may  be  applied  in  two  ways:  (1)  first  apply  BC's  at  x=-a, 
x=a  then  at  y=-b,  y=0-,  or  (2)  first  apply  BC's  at  y=-b,  y=0  then  at 
x=-a,  x=a. 

The  first  approach  may  be  interpreted  as  representing  the  Green's 

functions,  Gj,  interms  of  the  parallel  plate  waveguide  modes  propagating 
.  h 

in  the  ±  y  directions  and  reflecting  from  the  walls  at  y=-b,  y=0.  The 

second  approach  may  be  interpreted  as  representing  the  G5  in  terms 

h  . 

of  the  parallel  plate  waveguide  modes  propagating  in  the  ±  x  directions 
and  reflecting  from  the  walls  at  x=-a,  x=a.  For  reasons  of  convergence 


16 


(discussed  in  Chapter  III),  the  latter  approach  is  used  in  deriving  the 
cavity  Green's  functions.  The  resulting  equations  are 


-1  *  csc2Ya 

G-(x,y;x' ,y' )  =  b  ).  y  [cosY(2a-|x-x' |  )-cosy(x+x' )1 
^  1=1 


■  tiry' 

Sin  b  sin  5 


(2.12a) 


-1  “  csc2Ya 

Gh(x,y;x' ,y' )  =  “b  I  [cosY(2a-|x-x' I  )+cosy(x+x' )1 

£=0 


£ity  iny* 

cos  b  COS  b 


{2.12b) 


where 


-  I  bJ 


bk 

<  v 


'■se(s 


r  r 


>  Ti  ’^e(G^y^) 


and 


z 


to 


2 

1 


I  =  0 

a  0 


(2.13a) 

(2.13h) 


(2.13c) 


Note  that  the  other  representation  of  the  Green's  function  is  entirely 
equivalent  to  Equation  (2.12),  hut  is  simply  not  as  efficient  for  the 
geometry  and  method  of  solution  presented  here. 

Before  discussing  the  method  of  solution  to  the  integral  equations, 
a  note  should  be  made  concerning  the  aperture  field  edge  conditions 
inherent  in  a  geometry  such  as  that  of  Figure  2.3.  These  edge 
conditions  are  significant  since  a  solution  for  the  aperture  fields 
should  at  least  approximately  meet  them.  It  is  noted  that  the  edge 
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behavior  can  significantly  affect  the  rest  of  the  aperture  field 
distribution;  especially  for  grazing  and  near-grazing  angles  of 
incidence. 

Hurd  Clfil  considers  the  problem  of  field  edge  behavior  for  a 
conducting  wedge  in  the  presence  of  dielectric  wedges.  If  we  are 
interested  in  the  edge  behavior  of  the  aperture  electric  fields  for  the 
geometry  of  Figure  2.3,  then  the  geometry  of  Figure  2.4  should  be 
relevant  (assuming  a  non-vanishing  cavity).  For  this  geometry  Hurd 
derives  the  edge  conditions 


=  Cj  (2.14a) 

Ep  =  tp^"^  (2.14b) 


where 


(2.15) 


Thus,  for  finite  e^,  E^  vanishes  at  the  edge  (p=0)  and  E^  becomes 
infinite  at  the  edge. 


Figure  2.4.  Geometry  for  consideration  of  edge  conditions. 
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III.  SOLUTION  OF  INTEGRAL  EQUATIONS  BY  THE  MOMENT  METHOD 

A.  THE  MOMENT  METHOD 

The  moment  method  [1]  is  a  powerful  and  frequently  used  numerical 
technique  for  solving  integral  and  integro-differntial  equations.  It  is 
a  method  that  transforms  the  integral  equation  into  a  matrix  equation 
which  is  then  solved  by  matrix  inversion. 

The  integral  equations  (2.7)  and  (2.8)  are  of  the  form 

b 

/  G(x,x ' )  M(x ' )  dx '  =  F(x)  , 

a  (3.1) 

a  <  x,x'  <  b  , 

where  G  is  a  known  function  or  operator  which  constitutes  the  kernel  of 
the  integral  equation;  F  is  a  known  forcing  function,  and  M  is  the 
unknown  function  to  be  determined.  Note  that  for  the  problem 
considered,  M  is  an  equivalent  magnetic  current,  F  is  an  incident  field 
(or  its  derivative),  and  G  is  related  to  the  Green's  function.  The 
first  step  in  the  moment  method  solution  to  Equation  (3.1)  is  to  expand 
the  unknown  function  M(x')  in  a  series  of  basis  functions  (also  called 
expansion  functions)  with  unknown  coefficients.  The  basis  functions  may 
be  either  subsectional ,  in  which  each  function  spans  a  subsection  of  the 
domain  of  M(x'),  or  they  may  be  entire,  in  which  each  function  spans  the 
the  entire  domain  of  M(x').  Some  examples  of  subsectional  basis 
functions  are  pulses  and  piecewise-continuous  sinusoids.  Examples  of 
entire  basis  functions  include  Fourier  series  and  other  modal 
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expansions,  traveling  waves,  and  functions  derived  from  the  geometrical 
theory  of  diffractions  (GTD)  [6]  and  its  uniform  version  (UTD)  [14].  In 
this  work  we  expand  in  a  series  of  entire  basis  functions: 

N 

M(x')  D  H  (x')  ,  a  <  x'  <  b  (3.2) 

n=l  "  " 

where  Mp(x)  is  a  basis  function  and  is  the  corresponding  unknown 
coefficient.  Note  that  as  N  ■*•  •,  (3.2)  can  yield  a  formally  exact 
representation  of  M(x').  For  N  finite  but  large  the  representation  is 
approximate,  but  nearly-exact  results  are  possible.  Substituting  (3.2) 

into  (3.1)  and  assuming  that  (3.2)  is  a  uniformly  convergent  expansion, 

N  b 

I  0^  /  G(x,x')  M^(x’)  dx*  =  F(x)  .  (3.3) 

n=l  a 

As  it  stands.  Equation  (3.3)  is  a  single  equation  in  N  unknowns. 

To  convert  it  to  a  set  of  N  simultaneous  equations  in  N  unknowns,  a  set 
of  N  weighting  functions  is  chosen  and  applied  to  Equation  (3.3)  at  N 
locations  in  the  domain  of  the  weighting  function.  Thus  the  second  step 
in  the  moment  method  solution  is  the  selection  of  a  set  of  weighting 
functions  (also  called  testing  functions). 

Let  {w(x),  a  <  X  <  b|  be  a  set  of  N  weighting  functions.  Then,  we 
may  write  the  weighted  equations  as 

N  b  b 

I  On  /  /  G(x,x')  M^(x*)  W|^(x)  dx’  dx 
n  =  l  a  a 

b 

=  J  F(x)  W^{x)  dx  ,  m  =  1 . N.  (3.4) 

a 
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Examples  of  simple  and  effective  weighting  schemes  are  point -matching 
and  Galerkin's  method.  In  point -matching  the  weighting  functions  are 
simply  delta  functions, 

“m**'  "  ^  ‘  ’Sn  "  •  (3-5) 

In  Galerkin's  method  the  weighting  functions  are  identical  to  the  basis 
functions,  i.e., 

W^(x)  =  M^(x)  .  (3.6) 

Equations  (3.5)  and  (3.6)  separately  constitute  the  only  choices  of 
weighting  functions  considered  in  this  work. 

Equations  (3.4)  may  be  expressed  in  matrix  form  as 

['mnl  [Onl  '  <3.7) 

where 

b  b 

=  /  /  G(x,x')  M^(x')  W^(x)  dx’  dx  (3.8) 

a  a 

Vm  =  F(x)  W^(x)  dx  .  (3.9) 

a 

Note  that  is  the  so-called  impedance  matrix  and  [V^^]  is  the 

voltage  matrix.  Finally,  the  unknown  coefficients  are  determined  from 


The  accuracy  of  the  results  depends  on  the  choice  of  both  basis 
functions  and  weighting  functions.  The  basis  functions  should  be  chosen 
so  that  Equation  (3.2)  can  approximate  the  unknown  current  M(x')  fairly 
well.  The  edge  conditions  and  6(x,x‘)  may  place  restrictions  on  the 
useful  choice  of  basis  functions.  The  same  restrictions  may  apply  to 
the  weighting  functions,  fortunately  the  moment  method  is  a  variational 
procedure  [11  and  often  yields  good  results  for  a  wide  range  of  basis 
and  weighting  functions. 

The  moment  method  is  generally  considered  to  be  a  low-frequency 
technique  limited  at  present  to  objects  <  10  wavelengths  for  two 
dimensional  problems.  This  limit  is  a  computational  one  and  is  due  to 
the  large  number  of  unknowns  required  to  accurately  model  the  currents 
on  large,  complicated  objects.  For  subsectional  basis  functions  (such 
as  pulses)  typically  10-20  unknowns  per  wavelength  are  neessary. 
Appropriate  entire  basis  functions  can  reduce  this  by  a  factor  of  2  or 
more.  When  the  number  of  unknowns  exceeds  several  hundred,  the  time  and 
cost  required  to  invert  the  impedance  matrix  becomes  exhorbitant  and  the 
technique  is  no  longer  attractive. 

In  problems  where  one  can  identify  the  dominant  physical  mechanisms 
contributing  to  the  unknown  current,  a  hybrid  moment  method  [101, [111, 
can  be  a  powerful  tool.  Basis  functions  are  chosen  which  incorporate 
the  mathematical  form  of  the  known  physical  mechanisms.  More 
traditional  basis  functions  (e.g.  pulses)  may  be  used  in  locations 
where  the  physics  is  not  so  well  understood.  The  use  of  these  hybrid  or 
traveling  wave/RTD  basis  functions  can  yield  efficient  high-frequency 


22 


solutions  at  a  cost  of  increased  mathematical  complexity.  Many 
excellent  examples  of  this  hybrid  or  traveling  wave/GTD  approach  to  the 
moment  method  exist  in  the  literature  as  mentioned  earlier, 

B.  SOLUTION  FOR  RECTANGULAR  CAVITIES  OF  NARROW 
TO  MODERATE  WIDTH  AND  ARBITRARY  DEPTH 

In  this  section  we  present  a  solution  to  the  problem  depicted  in 
Figure  2.3  which  is  valid  for  either  loaded  or  unloaded  rectangular 
cavities  of  narrow  to  moderate  width  and  arbitrary  depth.  This  solution 
requires  approximately  four  unknowns  per  (free-space)  wavelength  across 
the  aperture,  and  is  therefore  limited  to  cavities  of  width  ~  25X  or 
less. 

Entire  basis  functions  are  used  with  a  Galerkin  weighting  scheme. 
The  basis  functions  are  chosen  to  be  the  transverse  portion  of  the 
modes,  both  propagating  and  evanescent,  of  a  parllel  plate  waveguide 
with  plates  at  x=-a  and  x=a.  This  choice  of  basis  functions  was  also 
made  by  Wang  [4]  in  one  part  of  his  solution  to  the  scattering  by  a  slot 
in  a  thick  conducting  screen. 

The  modes  supported  by  a  parallel  plate  waveguide  are  easily 
determined  from  the  2-D  scalar  wave  equation  and  the  boundary  conditions 
at  the  walls.  Expanding  the  unknowns  in  the  integral  Equations  (2.7) 
and  (2.8)  in  terms  of  the  transverse  portion  of  these  modes,  we  obtain 
for  the  hard  and  soft  polarizations,  respectively, 

3ufi  N  nn 

D  cos  ?!  (a-x')  (3.11) 

n=0  " 
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(3.12) 


I;  nw 

u.(x')  =  I  Dp  sin  ^  (a-x') 
n=l 

for  -a  <  x'  <  a 


Note  that  the  in  (3.11)  and  (3.12)  are  in  general  different. 

Applying  Galerkin's  method  we  determine  the  impedance  and  voltage 
matrices  from  Equations  (3.8)  and  (3.9).  For  the  case  of  hard 
polarization  the  results  are 

,  nir  m 

Kn  =  1  i  f’(x.x  )  cos  (a-x')  cos  (a-x)  dx'dx  (3.13) 


-a  -a 


m 


=  -2H  I  eJKxcos*' 


mn 

cos  23  (a-x)  dx 


(3.14) 


-a 


where  G(x,x')  equals  the  quantity  in  brackets  on  the  LHS  of  Equation 

(2.7) .  Similarly,  the  results  for  the  soft  polarization  are 

^mn  =  /  /  G(x,x')  sin  Ja  (a-x')  sin  Ja  (a-x)  dx'dx  (3.15) 

-a  -a 

=  2jksin<.'  Eg  j%jkxcos^'  ^  ^ 

•  3 

Where  G(x,x')  equals  the  quantity  in  brackets  on  the  LHS  of  Equation 

(2.8) .  Note  that  the  impedance  matrix  is  symmetric. 

The  integrals  in  the  expressions  for  are  easily  evaluated; 
however,  the  double  integrals  in  the  expressions  for  Z^n  are 
considerably  more  challenging.  The  mathematical  details  of  this 
evaluation  are  given  in  Appendix  C.  The  double  integrals  of  the  Hankel 
function  terms  in  the  total  Green's  functions  are  reduced  to  single 
integrals  by  a  rotation  of  coordinates  (the  so-called  Popovich 
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transformation).  The  single  integrals  are  evaluated  over  the  initial 
part  by  numerical  integration  and  the  remaining  part  is  evaluated 
analytically  by  using  the  large-argument  form  of  the  Hankel  function. 

The  double  integrals  of  the  cavity  Green's  functions,  given  by 
Equations  (2.12)  and  (2.13),  can  be  evaluated  in  closed  form,  leaving 
infinite  sums.  Specifically  for  this  form  of  the  cavity  Green's 
functions,  portions  of  these  infinite  sums  can  be  summed  into  closed 
form.  For  the  more  conventional  choice  of  Green's  functions,  as  given 
by  Equations  (B.22)  and  (B.23),  such  closed-form  terms  were  not 
obtained.  This  fact  accounts  for  the  significantly  faster  convergence 
for  the  chosen  cavity  Green's  functions.  Only  a  few  evanescent  terms 
are  required  to  achieve  numerical  convergence  of  the  remaining  portion 
of  the  sum. 

We  point  out  that  the  current  representation  given  by  Equation 
(3.12)  satisfies  the  edge  condition  (2.14a)  in  the  sense  that  (3.12) 
goes  to  zero  for  [x'l+a.  Actually  each  basis  function  in  (3.12) 
vanishes  at  the  edges  of  the  aperture,  a  fact  which  greatly  simplifies 
the  evaluation  of  Z^n-  For  the  hard  polarization  the  current 
representation  given  by  Equation  (3.11)  will  not  go  to  infinity  at  the 
aperture  edges  as  required  by  (2.14b),  assuming  a  finite  value  of  N. 
However,  (3.11)  will  approximate  the  edge  condition  well  enough  to  yield 
accurate  results  even  with  only  four  unknowns  per  wavelength. 

Finally,  the  far-zone  scattered  fields  in  region  I  (y>0)  are  found 
from  Equations  (2.9),  (3.11)  and  (3.12)  to  be 
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C.  SOLUTION  FOR  RECTANGULAR  CAVITIES 
OF  LARGE  WIDTH  AND  SMALL  DEPTH 

In  this  section  we  present  a  solution  to  the  problem  depicted  in 
Figure  2.3  which  is  valid  for  loaded  rectangular  cavities  of  large  width 
and  small  depth.  The  number  of  unknowns  in  this  solution  is  independent 
of  the  width  of  the  aperture,  but  can  depend  on  the  depth  of  the  cavity. 
By  small  depth  it  is  meant  that  the  dielectric  portion  of  Figure  2.3 
will  support  at  most  one  surface  wave  mode. 

In  the  previous  section,  a  solution  involving  a  current  expansion 
in  parallel  plate  waveguide  modes  was  presented.  This  solution  will 
also  work  for  the  present  case  of  wide,  shallow  cavities,  but  it  will 
require  many  unknowns  (3-4  per  wavelength  across  the  aperture).  By 
considering  the  physics  of  the  wide,  shallow  cavity  we  can  arrive  at  a 
more  efficient  set  of  basis  functions. 

Consider  the  grounded  dielectric  slab  shown  in  Figure  3.1.  This 
structure  can  support  surface  waves,  incident  and  reflected  geometrical 
optics  fields,  plus  other  types  of  waves  which  we  are  not  concerned  with 
here.  For  the  case  of  plane  wave  incidence  the  total  field  in  the 
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region  above  the  slab  is  simply  the  incident  plus  reflected  geometrical 
optics  terms  (since  a  plane  wave  cannot  excite  surface  waves  in  an 
infinite  dielectric  slab): 


Us  (x.y)  =  {  'o  I  •  [ejk(yrt)sfn*'  ^  ^-jk(y*b)sin*'j 

f;  h 

.  gjkxCOS(J)' 

where  the  reflection  coefficient  is  [191 
j  sin(|)'  +  h  •  tan  (kbh) 

Y  = - 

^  j  Gp  Up  sin(ji'  -  h  •  tan  (kbh) 
j  sin^i'  -  -  .ot  (kbh) 

r  -  - -  -  ■  - 

s  j  sin4i'  -  n  •  cot  (kbh) 
and 

h  -  /  Cp  Up  -  sin^  i|)' 
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GROUND  PLANE 


Figure  3.1.  Scattering  from  an  infinite  dielectric  slab  on 
plane. 
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Note  that  the  reflection  coefficient  has  a  magnitude  of  unity,  as  we 
would  expect  for  an  incident  plane  wave. 

If  we  truncate  the  slab  of  Figure  3.1,  then  the  diffration  of  an 
incident  plane  wave  at  the  ends  of  the  slab  will  excite  surface  waves. 
The  z-components  of  hard  and  soft  normalized  surface  waves  are  [161 
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where  'hard'  and  'soft'  are  as  defined  by  Equation  (2,1).  The  surface 
wave  propagation  constant,  3,  is  determined  from 

Oj  tan  a^b  *  (3.23a) 


-a,  cot  a.b  = 

1  1  r 


(3.23b) 


for  the  hard  and  soft  polarizations,  respectively,  where 

a2  =  b2  .  k2 

“l^  = 


(3.24a) 

(3.24b) 


Substituting  Equations  (3.24)  into  (3.23)  yields  a  transcendental 
equation  in  B,  which  is  then  solved  by  the  Newton-Raphson  method.  From 
Equations  (3.23a ,b)  it  may  be  shown  that  the  first  surface  wave  mode  for 
the  hard  polarization  has  no  cutoff  frequency.  This  is  not  true  for  the 


soft  polarization.  Thus  even  a  very  thin  grounded  dielectric  slab  will 
support  a  surface  wave  only  for  a  hard-polarized  incident  field. 

We  now  return  to  the  geometry  of  Figure  2.3  and  consider  a  suitable 
choice  of  basis  functions.  From  the  above  discussion  it  seems 
reasonable  to  assume  that  a  geometrical  optics  type  of  field  (also 
called  a  forced  wave  inside  the  dielectric),  such  as  Equation  (3.18), 
will  make  a  significant  contribution  to  the  aperture  field.  In 
addition,  if  we  assume  the  dominant  surface  wave  mode  is  present,  then 
terms  such  as  (3.21)  and  (3.22)  are  also  important.  Thus  far,  our 
expansion  for  the  aperture  electric  fields  is 

=  Dj  +  D3  e-^^^  ;  |x|<a  (3.25) 

where  8  depends  on  polarization.  The  coefficients  ^2  ^3  include  the 

multiple  reflections  of  the  surface  waves  from  the  end  walls  at  x=-a  and 
x=a.  Richmond  [12]  employed  a  very  similar  expansion  in  his  moment 
method  solution  to  scattering  by  a  thin  dielectric  strip  (E^ 
polarization  only). 

For  this  problem  considered  here  the  current  representation  (3.25) 
does  not  properly  account  for  the  aperture  physics  near  the  end  walls. 
The  transition  region  exists  between  the  edges  of  the  aperture  and  the 
region  where  the  surface  waves  become  established.  This  transition 
region  typically  extends  about  one  wavelength  into  the  aperture  from  the 
edges,  except  for  near  grazing  angles  of  incidence  where  it  could 
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possibly  extend  over  the  entire  aperture  region  if  the  aperture  is  only 
moderately  wide. 

To  account  for  the  physics  of  the  transition  regions  additional 
basis  functions  were  chosen  which  were  entire  but  of  importance  only 
near  the  edges.  To  keep  the  mathematical  complexity  to  a  minimum,  a  set 
of  exponential  functions  which  decay  away  from  the  edges  were  chosen. 
Thus  the  current  expansion  now  becomes 
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(3.26) 


where  the  q^  are  empirically  determined  complex  nuidbers  with  a  negative 
imaginary  component.  A  discussion  of  suitable  values  for  the  q^-  is 
given  in  the  next  chapter;  typically  there  are  several  transition  terms. 
Note  that  the  surface  waves  and  the  transition  terms  in  Equation  (3.26) 
are  written  in  a  form  that  indicates  that  these  terms  emanate  from  the 
edges  of  the  aperture,  whereas  the  forced  wave  term  is  not  written  in 
such  a  form. 

Recall  that  U5(x)  must  satisfy  the  edge  condition  (2.14a).  Thus, 
for  the  soft  polarization.  Equation  (3.26)  must  be  set  equal  to  zero  at 
x=-a  and  x=a.  iJpon  doing  this  we  may  then  solve  for  two  of  the  N 
unknown  coefficients  in  terms  of  the  remaining  coefficients,  thus 
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reducing  the  nurt)er  of  unknowns  from  N  to  N-2.  The  resulting  form  for  a 
given  basis  function  for  the  soft  polarization  becomes 


-jO(a-x)  ^  ^  ^  _  ,3  2,, 


Mr,(x)  =  a„  e 


where  hp  and  fp  are  the  appropriate  values  of  kcoS(p,  8,  and  shown  in 
Equation  (3.26),  and  0  is  an  empirically  determined  complex  number  with 
a  negative  imaginary  component  (like  the  q^-),  and  where 
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A  value  of  0  is  chosen  so  that  the  basis  function  given  by  Equation 
(3.27)  goes  to  zero  rapidly  near  the  edges  of  the  aperture;  a  discussion 
of  actual  numerical  values  is  deferred  to  the  next  chapter.  The 
corresponding  lorm  of  a  given  basis  function  for  the  hard  polarization 
i  s 
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Apnlying  Ralerkin's  method  we  determine  the  impedance  and  voltage 
mat"ices  from  Equations  (3.R)  and  (3.9).  For  the  hard  polarization  the 
results  are 
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where  G(x,x')  equals  the  quantity  in  brackets  on  the  LHS  of  Equation 
(2.7).  Similarly  the  results  for  the  soft  polarization  are 


"3  “3 
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+  dx  (3.33) 

where  G(x,x')  equals  the  quantity  in  brackets  on  the  LHS  of  Equation 
(2.8).  Note  that  the  impedance  matrix  is  symmetric. 

Again  the  expressions  for  V^,  are  easily  evaluated.  Comparing  the 
above  expressions  for  the  impedance  matrix  elements  with  those  of  the 
waveguide  mode  current  expansion,  given  by  Equations  (3.13)  and  (3.15), 
we  see  that  in  both  cases  the  basis  and  weighting  functions  may  be 
expressed  as  exponentials.  Thus,  the  evaluation  of  Equations  (3.30)  and 
(3.32)  is  not  fundamentally  different  from  that  of  Equations  (3.13)  and 
(3.15). 

Finally  the  far-zone  scattered  fields  in  region  I  (y>0)  are  found 
from  Equations  (2.9),  (3.27)  and  (3.29)  to  be 
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for  the  soft  polarization,  and 
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for  the  hard  polarization. 
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IV.  NUMERICAL  RESULTS  AND  DISCUSSION 


A.  CAVITIES  OF  ARBITRARY  DEPTH  AND  NARROW  TO  MODERATE  WIDTH 
RECESSED  IN  A  GROUND  PLANE 

In  this  section  we  consider  the  EM  scattering  from  loaded  or 
unloaded  cavities  of  narrow  to  moderate  width  and  recessed  an  arbitrary 
distance  in  a  ground  plane.  By  "narrow  to  moderate  width"  it  is  meant 
that  the  aperture  width  is  <  10  X,  where  X  is  the  free-space  wavelength. 

Before  discussing  the  numerical  results  in  this  section,  it  is 
important  to  introduce  an  explanation  of  the  titles,  labels  and  variable 
names  that  accompany  each  of  the  computer  plots  pertaining  to  the 
numerical  calculations  presented  below. 

TITLES  OF  PLOTS 

PMH  =  pulses/point  matching  solution  for  hard  polarization 

NOTCHH,  NOTCHS  =  hybrid  basis  solution  for  hard,  soft 
polarization 

(ignore  the  .TEST  appearing  in  some  titles) 

COORDINATE  LABELS 

BACKSCATTER  =  monostatic  backscattered  field  or  H^) 

magnitude  multiplied  by  the  factor  /nTT,  where  k  is 
the  free-space  wavenumber. 

MAGNITUDE,  PHASE  =  equivalent  current  magnitude  and  phase 

PATTERN  MAGNITUDE  *  bistatic  scattered  field  (E^  or  H^) 
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VARIABLES  COMMON  TO  ALL  SOLUTIONS 

WIDTH  =  aperture  width,  X  (free-space) 

DEPTH  =  cavity  depth,  X 

ERR  =  Re  (cp)  in  the  dielectric  loading  of  the  cavity 

ERI  =  Im  (cp)  in  the  dielectric  loading  of  the  cavity 

lIRR  =  Re  (up)  ifi  the  dielectric  loading  of  the  cavity 

UR  I  =  Im  (up)  in  the  dielectric  loading  of  the  cavity 

PHP  =((>',  0°  <(})'  <  90°  (it>'=90°  =>  normal  incidence.) 

VARIABLES  SPECIFIC  TO  A  SOLUTION 

PMH  =  MM  solution  using  pulse  expansion  and  point  matching 
in  the  aperture  only  for  the  hard  case 

NMODE  =  Number  of  modes  in  the  cavity  Green's  function 

A 

(modes  propagating  in  ±  y  direction) 

PW  =  Pulse  width,  X 

N1,N2  =  Nunber  of  pulses  counting  from  x=“a,  x=+a. 


NOTCHH.NOTCHS  =  Hard  and  soft  case  MM  solutions  using  modal 
expansion  in  waveguide  aperture  respectively 

NZ1,NZ  =  Nurber  of  unknown  modal  basis  functions  for  hard, 
soft  polarization 

NMODE  =  Nunber  of  terms  in  reduced  form  of  cavity  Green's 
function  (modes  propagating  in  ±  x) 

DX  =  2  *  integration  step  of  numerical  portion  of 
integration,  X 
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EUAVH3,EUAVS  *  Hard  and  soft  case  MM  solution  using  a  hybrid 

traveling  wave  expansion  in  the  aperture,  respectively. 

NMOOE  =  Same  as  for  NOTCHH.NOTCHS 

OX  =  Same  as  for  NOTCHH.NOTCHS 

NFW  =  Number  of  forced  wave  basis  functions 

NSW  =  1/2  number  of  surface  wave  basis  functions 

NEW  =  1/2  number  of  transition  term  basis  functions 

We  first  consider  results  for  the  hard  polarization.  Figure  4.1 
compares  a  pulse/point -matching  solution  for  the  monostatic  scattered 
magnetic  field  (H^)  from  a  loaded  recessed  cavity  of  width=lx, 
depth=0,lx  to  the  fialerkin  solution  employing  entire  waveguide  mode 
basis  functions  (Chapter  III,  Section  B).  The  pulses/point-matching 
solution  employs  20  pulses  of  width=0,05X  and  10  modes  in  the  cavity 
Green's  function  (using  the  form  given  by  Equation  (R.22b);  the  CPU  time 
on  a  PRIME  computer  was  49  seconds.  The  modal  basis  solution  employs  5 
waveguide  modes,  no  modes  in  the  expression  for  Zl  of  (C.12)  (the 
closed-form  portion  of  Zi_  proved  to  be  sufficient),  and  a  numerical 
integration  step  of  0.05X;  the  CPU  time  was  15  seconds.  Note  that  both 
solutions  were  verified  to  have  converged.  Similar  results  are  given  in 
Figure  4.2  for  a  depth=10.1X.  Again,  there  is  good  agreement  between 
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Figure  4.2.  Backscattering  by  a  deep  notch  cavity  in  a  ground  plane 
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the  solutions.  Note  that  even  for  the  deep  cavity,  no  modes  were 
required  for  the  computation  of  Zl* 

Now  consider  the  soft  polarization.  Figures  4.3  and  4.4  show  the 
monostatic  scattered  electric  field  (E^)  computed  from  the  modal  basis 
method  for  cavities  of  depth  =  O.IA  and  10. lA,  respectively.  This 
solution  employs  4  waveguide  modes,  no  modes  in  the  expression  for 
(again,  the  closed-form  portion  of  Z|_  was  sufficient),  and  an 
integration  step  of  .05A.  Bistatic  scattered  field  equivalent  currents 
and  far-field  patterns  computed  with  the  modal  basis  solution  were 
compared  with  a  moment  method  solution  [4]  employing  piecewise 
continuous  sinusoids  and  Galerkin's  method.  Good  agreement  was  observed 
for  the  two  solutions  for  normal  incidence.  There  appears  to  be  some 
question  as  to  the  use  of  the  computer  codes  in  [4]  for  grazing  angles; 
hence,  this  comparison  has  not  been  made  at  this  time. 

B.  LOADED  SHALLOU  CAVITIES  OF  LARGE  WIDTH 
RECESSED  IN  A  GROUND  PLANE 

In  this  section  we  consider  scattering  from  loaded  shallow  cavities 
of  large  width  (>  5  A)  recessed  in  a  ground  plane.  By  "shallow"  it  is 
meant  that  at  most  one  surface  wave  is  supported  in  the  dielectric. 

This  condition  places  constraints  on  the  depth  of  the  cavity  and  the 
electrical  properties  of  the  dielectric  medium  according  to  Equations 
(3.23)  and  (3.24).  Table  4.1  shows  surface  wave  propragation  constants 
as  a  function  of  dielectric  thickness  (i.e.,  cavity  depth)  for  a 
dielectric  with  ep2.5,  yp=i.o.  A  dash  (-)  indicates  that  a  surface  wave 
is  not  supported. 
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Figure  4. 
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Figure  4.4.  Backscattering  by  a  deep  notch  cavity  in  a  ground  plane 


TABLE  4.1 


HARD  AND  SOFT  SURFACE  HAVE  PROPAGATION  CONSTANTS 
FOR  A  GROUNDED  DIELECTRIC  WITH  e 
AS  A  FUNCTION  OF  DIELECTRIC 
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6.789 

— 

.15 
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— 

.20 
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6.288 
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8.599 

6.630 

Figure  4.5  compares  a  pulses/point-matching  solution  to  the  modal 
basis  solution  for  the  monostatic  scattered  magnetic  field  from  a  loaded 
recessed  cavity  of  width=5A,  depth=.05A,  ep=2.5.  The 
pulses/point -matching  solution  employs  50  pulses  of  width=.10A  and  40 
modes  in  the  cavity  Green's  function;  the  CPU  time  was  9  minutes,  47 
seconds.  The  modal  basis  solution  employs  21  waveguide  modes,  only  the 
closed-form  term  in  Zl,  and  an  integration  step  of  .05A;  the  CPU  time 
was  1  minute,  10  seconds.  The  modal  basis  solution  was  observed  to 
converge.  Due  to  lack  of  disk  space  we  were  unable  to  verify  the 
convergence  of  the  point  matching  solution. 


r=2.5,  vir=l»0 


r- vip- 

THICKNESS 
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Also  shown  in  Figure  4.5,  at  ten  degree  intervals,  is  the  solution 
obtained  using  the  basis  functions  discussed  in  Section  C  of  Chapter 
III.  This  solution  employs  a  forced  wave,  two  surface  waves,  and  8 
transition  terms  for  a  total  of  11  basis  functions,  and  is  valid  for  any 
cavity  width  larger  than  a  few  wavelengths.  The  surface  wave 
propagation  constant  for  this  geometry  is  given  in  Table  4.1.  The 
transition  term  "propagation"  constants  were  chosen  to  be  of  the  form 
qi  =  k  -  js^  .  (4.1) 

where  the  are  defined  in  Equation  (3.26),  k  is  the  free-space 
wavenumber,  and  the  s,-  are  empirically  determined  attenuation  constants. 
It  was  found  from  numerical  experimentation  that  allowing  the  Sj  to  take 
on  the  four  integer  values 

»  5,6. 7,8  .  (4.2) 

resulted  in  accurate  monostatic  scattered  fields  for  shallow  cavities. 
Actually,  the  nonostatic  scattered  field  is  somewhat  insensitive  to  the 
specific  q-j  chosen,  as  long  as  the  corresponding  basis  functions  are 
insignificant  a  wavelength  or  so  from  the  aperture  edge.  The  CPU  time 
for  the  points  shown  in  Figure  4.5  was  1  minute,  16  seconds,  which  is 
similar  to  the  CPU  required  for  the  modal  basis  solution.  Even  though 
fewer  unknowns  are  required  in  the  hybrid  basis  solution  (11  versus  21 
in  the  modal  basis  solution),  one  of  the  basis  functions,  the  forced 
wave,  is  a  function  of  the  angle  of  incidence,  as  can  be  seen  from 
Equation  (3.26).  Thus,  a  row  of  the  impedance  matrix  must  be  recomputed 
and  the  resulting  matrix  inverted  for  each  new  angle  of  incidence.  In 
the  modal  basis  solution  the  impedance  matrix  is  computed  and  inverted 
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only  once.  The  advantage  of  the  hybrid  basis  solution  is,  of  course, 
that  only  11  unknowns  are  required  regardless  of  cavity  width. 

Figures  4.6  and  4.7  compare  the  hybrid  and  modal  basis  solutions 
for  the  equivalent  magnetic  aperture  current  and  bistatic  scattered 
field  pattern  for  a  plane  wave  normally  incident  on  the  geometry  of 
Figure  4.5,  except  that  the  cavity  depth=0.10x.  The  required  CPU  for 
the  modal  basis  solution  was  54  seconds;  versus  28  seconds  for  the 
hybrid  basis.  Figures  4.8  and  4.9  make  the  same  comparison  for  grazing 
incidence  (!(i'=0°).  Note  that  excellent  agreement  between  the  two 
solutions  even  for  grazing  incidence. 

Figures  4.10  to  4.13  compare  the  monostatic  scattered  field 
obtained  from  the  hybrid  basis  solution  at  ten  degree  intervals  to  the 
modal  basis  solution  for  cavity  depths  of  .lOX,  .15X,  .20X,  and  .25  X. 
Good  agreement  is  obtained  everywhere  except  for  the  grazing  incidence 
region  for  a  depth  of  .25X,  the  largest  depth  which  was  tried.  Such 
variations  might  yield  good  results  for  a  particular  geometry,  but  would 
not  always  apply  to  other  geometries.  Thus,  a  consistent  method  for 
choosing  exponential -type  transition  terms  was  not  found  for  the  deeper 
cavities,  i.e.,  for  depth=.20X  and  .25X.  In  short,  an  accurate  and  very 
efficient  hybrid  solution  for  the  hard  polarization  was  found  for 
cavities  with  depths  less  than  .15X  and  of  arbitrary  width. 

Consider  now  the  soft  polarization  results  of  plane  wave  scattering 
from  the  same  geometry  considered  in  the  above  hard  polarization 
analysis.  From  Table  4.1  one  sees  that  surface  waves  do  not  exist  in 
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Figure  4.6.  Current  distribution  on  the  aperture  of  a  wide  and  shallow 
cavity  for  normal  incidence. 
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Figure  4.8.  Current  distribution  in  the  aperture  of  a  wide  and  shallow 
cavity  for  grazing  incidence. 
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the  thin  dielectrics  for  the  material  properties  ep=2.5,  iip=1.0. 

However,  a  transition  surface  wave  type  of  field  does  exist  under  cutoff 
conditions.  From  plots  of  the  equivalent  currents  in  the  aperture  it  is 
readily  apparent  that  this  field  has  the  propagation  constant  of  free 
space  and  becomes  more  significant  as  cutoff  of  the  dominant  surface 
wave  is  approached.  For  a  depth  of  .15X  and  .lOX  it  is  less  obvious. 

For  this  reason  a  "surface  wave"  with  propagation  constant=k  will  be 
employed  in  the  hybrid  basis  solution  for  the  recessed  cavity  of 
depth=.15x,  even  though  the  dominant  surface  wave  is  actually  cutoff. 

Figure  4.14  compares  the  hybrid  basis  solution  to  the  modal  basis 
solution  for  the  monostatic  scattered  electric  field  from  a  loaded 
recessed  cavity  of  width=5x  and  depth=.05x.  The  modal  basis  solution 
employs  20  modes;  the  CPU  time  was  1  minute  7  seconds.  Because  of  the 
fundamental  similarity  between  the  modal  basis  solutions  for  the  hard 
and  soft  polarizations,  the  CPU  times  are  nearly  identical  for  a  given 
number  of  unknowns.  The  hybrid  basis  solution  employs  only  a  forced 
wave.  However,  recall  that  additional  exponential  terms  are  employed  so 
that  each  basis  function  will  meet  the  edge  condition.  Thus,  for  a 
basis  function  of  the  form  given  by  Equation  (3.27),  a  value  of  0  must 
be  chosen  so  that  the  basis  function  will  go  to  zero  in  an  appropriate 
manner  as  (xl+a.  A  suitable  form  for  0  is 

0  =  k  -  js  .  (4.3) 

where  k  is  the  free-space  wavenumber,  and  s  is  determined  empirically. 
The  value  of  s  used  in  Figure  4.14  and  the  following  figures  is 

S  =  6  (4.4) 
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However,  nunerical  experimentation  has  shown  that  the  hybrid  basis 
solution  is  quite  insensitive  to  specific  value  of  0.  Any  choice  of  0 
which  takes  the  basis  function  to  zero  within  a  fraction  of  wavelength 
of  the  edge  is  suitable.  Note  that  the  exponential  transition  terms 
used  in  the  hybrid  basis  solution  for  the  hard  polarization  were  not 
found  to  be  useful  for  the  soft  polarization,  except  as  a  means  of 
enforcing  the  edge  condition. 

Because  there  is  only  one  unknown  in  the  hybrid  basis  solution 
shown  in  Figure  4.14,  the  CPU  time  is  only  20  seconds. 

The  agreement  between  the  hybrid  basis  and  modal  basis  solutions  in 
Figure  4.14  is  very  good  except  near  grazing  incidence.  This  disparity 
is  more  clearly  illustrated  in  the  next  set  of  figures.  Figures  4.15  to 
4.18  compare  the  equivalent  aperture  currents  and  bistatic  scattered 
field  patterns  of  the  two  solutions  for  angles  of  incidence  and 

(|i’=in°  and  a  cavity  of  depth  =  .inx.  Again,  only  the  forced  wave  is 
employed  in  the  hybrid  basis  solution,  and  20  modes  are  employed  in  the 
modal  basis  solution;  the  CPU  times  are  8  seconds  and  50  seconds 
respectively.  The  agreement  for  41 '=90°  is  very  good,  although  from 
Figure  4.15  one  can  see  that  a  weak  field  is  present  in  addition  to  the 
forced  wave  portion  of  the  total  field.  This  additional  field  is  more 
significant  for  *'=10°,  as  can  be  readily  seen  in  Figure  4.17.  The 
corresponding  bistatic  scattered  field  patterns,  shown  in  Figure  4.18, 
display  good  agreement  except  for  the  region  of  grazing  incidence, 
*'<30°.  Attempts  to  improve  the  solution  in  this  region  by  the  addition 
of  a  surface  wave  with  propagation  constant  k  and/or  the  inclusion  of 
exponential  transition  terms  were  unsuccessful. 
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Figure  4.15.  Current  distribution  in  the  aperture  of  a  wide  and  shallow 
cavity  for  normal  incidence. 
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Figures  4.19  to  4.22  show  the  monostatic  scattered  field  results 
for  cavities  of  depths  .lOX,  .15A,  .2nx  and  .25  X.  Inaccuracies  in  the 
grazing  incidence  region  are  plainly  visible  in  Figures  4.19  and  4.20. 
Inclusion  of  a  surface  wave  type  of  field  in  the  hybrid  solution  for 
b=.15X  improved  the  solution  which  employed  a  forced  wave  only, 
especially  for  normal  incidence.  However,  the  grazing  incidence 
solution  was  only  slightly  improved.  The  hybrid  solutions  for  b=.2nx 
and  .25X  employ  a  forced  wave  and  true  surface  waves,  since  the  dominant 
surface  wave  mode  is  no  longer  cutoff.  Figures  4.21  and  4.22  show  the 
poor  agreement  between  the  hybrid  and  modal  basis  solutions. 

Fxamination  of  the  bistatic  scattered  field  patterns  for  these 
geometries  shows  fairly  good  agreement  for  normal  incidence,  but  very 
poor  agreement  for  grazing  incidence.  Attempts  to  improve  the  hybrid 
solution  by  the  use  of  transition  terms  occasionally  yielded  better 
results.  However,  as  was  true  for  the  hard  polarization,  a  consistently 
accurate  set  of  basis  functions  was  not  found  for  the  deeper  cavities. 

In  short,  a  fairly  accurate  and  very  efficient  hybrid  solution  for  the 
soft  polarization  was  found  for  cavities  with  depths  <  .lOX  and  of 
arbitrary  width. 

An  alternative  hybrid  scheme  for  both  the  hard  and  soft 
polarizations  would  be  to  use  the  proper  transition  functions  associated 
with  the  phenomenon  of  edge  diffraction  within  the  surface  wave 
transition  region.  Such  functions,  involving  Fresnel  integrals  and 
terms  of  the  type  could  be  determined  from  asymptotic  study  of 

an  appropriate  canonical  problem.  However,  use  of  these  transition 
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functions  in  a  hybrid  moment  method  solution  would  greatly  increase  the 
CPU  required  to  compute  the  assuming  Galerkin's  method  is  employed. 
Use  of  point-matching  may  alleviate  this  problem. 

Next,  we  consider  the  effects  of  loading  on  the  scattering  patterns 
of  the  wide  shallow  cavities.  Comparisons  will  be  made  using  only 
results  from  the  modal  basis  solution.  Considering  first  the  hard 
polarization.  Figures  4.23  th'^ough  4.27  show  the  backscattered  field 
magnitude  for  a  cavity  of  width=5X  and  depth=.ix  and  for  various 
dielectric  loading;  in  all  cases  Up=l.  These  figures  should  be  compared 
with  Figure  4.10  for  which  €^=2.5  and  for  which  a  significant  surface 
wave  is  present.  The  backscatter  of  an  unloaded  cavity,  shown  in  Figure 
4.  .3,  shows  significantly  less  backscatter  along  grazing  angles  due  to 
the  lack  of  surface  wave.  A  lossy  dielectric,  ep=2.5-j.25,  is  shown  in 
Figure  4.24.  Again  the  most  significant  differences  compared  to  Figure 
4.10  occur  for  grazing  angles.  Figure  4.25  is  for  a  dielectric  constant 
Ep-S  and  displays  no  significant  differences  in  backscatter.  However, 

F  gures  4.26  and  4.27  are  for  larger  dielectric  constants,  ep=10  and 
ep=15,  respectively,  and  show  a  pronounced  decrease  in  backscatter  along 
g: izing  angles  as  compared  w’th  Figure  4.10.  This  effect  is  due  to  the 
destructive  interference  of  surface  wave  reflections  and  will, 
therefore,  vary  with  the  cavity  width. 
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Figure  4.23.  Backscattering  from  a  wide  and  shallow  cavity. 
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Figure  4.24.  Backscattering  fron  a  wide  and  shallow  cavity. 
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Figure  4.25.  Rackscattering  from  a  wide  and  shallow  cavity, 
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Figure  4.26.  Rackscattering  from  a  wide  and  shallow  cavity 
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Figure  4,27.  Backscatteri ng  from  a  wide  and  shallow  cavity. 


Figure  4.28  shows  the  backscatter  from  a  cavity  of  width=9.5\  and 
depth=.25X  with  a  dielectric  constant  ep=2.5.  This  geometry  is  easily 
shown  to  be  surface  wave  resonant  (3^=8. 599  from  Table  4.1);  however, 
this  is  the  case  is  readily  apparent  from  the  plot.  Figure  4.29 
dramatically  displays  the  effect  of  adding  loss  to  the  dielectric, 
ep=2.5-j.25,  thus  removing  the  resonance  by  both  altering  the  resonant 
frequency  and  damping  the  surface  waves. 

Consider  now  the  effects  of  loading  on  the  backscatter  patterns  for 
the  soft  polarization.  Figures  4.30  through  4.34  are  for  a  cavity  of 
width=5X  and  depth=.lx  and  for  various  types  of  dielectric  loading. 

These  should  be  compared  with  Figure  4.19  for  which  ep=2.5  and  for  which 
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Figure  4.?8.  Backscattering  from  a  wide  and  shallow  cavity. 


NflirilM 

«fTS 


WIPTM 

nii'Tii 

f  !■  ( 

1IMR 
llh  1 
I'llr 

■  1 
•WIPf 


nr  oK’i  I 


Figure  4.29.  Backscattering  from  a  wide  and  shallow  cavity. 
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a  significant  surface  wave  is  not  present.  Figures  4.30  and  4.31  show 
the  backscatter  of  an  unloaded  cavity  and  of  a  cavity  containing  a  lossy 
dielectric  with  ep=2.5-j.25,  respectively.  Neither  of  these  show  much 
difference  when  compared  with  Figure  4.19  simply  because  of  the  lack  of 
surface  wave  in  each  case.  The  same  rationale  applies  to  Figure  4.32 
for  which  ep=5,  even  though  the  surface  wave  is  only  slightly  below 
cutoff.  For  Figure  4.33,  for  which  ep=10,  the  surface  wave  mode  (or 
modes)  are  above  cutoff  as  is  readily  apparent  from  the  backscatter 
pattern.  Figure  4.34,  for  which  ep=15,  shows  a  decrease  in  backscatter 
compared  with  Figure  4.33  for  near-grazing  angles,  although  not  nearly 
as  pronounced  as  was  observed  for  the  hard  polarization. 

Figure  4.35  shows  the  backscatter  from  an  unloaded  cavity  of 
width=5X  and  depth=.25X.  Comparing  this  with  Figure  4.22,  for  which 
ep=2.5  and  for  which  a  surface  wave  is  present,  the  significance  of  the 
surface  wave  at  grazing  and  near-grazing  angles  is  readily  apparent. 
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Figure  4.30.  Backscatteri ng  fron  a  wide  and  shallow  cavity 
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Figure  4.31.  Backscatteri  ng  froi^  a  wide  and  shallow  cavity 
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NOTCHS 


Figure  4.32.  Rackscattering  from  a  wide  and  shallow  cavity. 


•  tmfi  9  oo 
|M  rxii  O  Ml 
»Mit  in  (Ml 
rH(  o  (Ml 
iihfl  1  on 
iM.i  n  'Ml 

l-llf  90  (Ml 

«»/  ;;o 

Mwrioc  O 

III  O  OS 


NDFCMS 


Figure  4.33.  Backscatteri ng  from  a  wide  and  shallow  cavity. 
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Figure  4.34.  Backscatttri ng  from  a  wide  and  shallow  cavity. 


Figure  4.35.  Backscattering  from  a  wide  and  shallow  cavity. 
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V.  CONCLUSIONS 


Two  efi'icient  moment  method  solutions  have  been  developed  for  the 
scattering  of  an  electromagnetic  plane  wave  from  a  2-0  rectangular 
cavity  recessed  in  a  perfectly  conducting  ground  plane  and  loaded  with  a 
homogeneous  lossy  or  lossless  dielectric.  Integral  equations  for  the 
equivalent  magnetic  current  in  the  cavity  aperture  are  solved  by 
Galerkin's  method  for  the  hard  and  soft  polarizations  independently. 

One  of  the  solutions  termed  a  modal  basis  solution,  expands  the 
equivalent  magnetic  current  in  a  set  of  parallel  plate  waveguide  modes 
which  constitute  entire  domain  basis  functions  over  the  whole  aperture. 
This  solution  requires  3-4  unknowns  per  aperture  wavelength,  versus 
10-20  unknowns  per  wavelength  when  using  pulses  or  piecewise  sinusoids. 
In  addition,  a  portion  of  the  cavity  Green's  function  contribution  to 
the  impedance  matrix  elements  was  summed  into  closed  form  resulting  in  a 
highly  convergent  representation.  When  compared  with  a 
pulses/point-matching  solution  for  the  hard  polarization,  the  modal 
basis  solution  was  found  to  be  several  times  as  efficient  without  loss 
of  accuracy.  At  present,  this  solution  is  limited  to  geometries  with 
aperture  widths  <  25X  (requiring  75-100  unknowns)  as  a  matter  of 
computational  efficiency;  the  required  CPU  is  independent  of  cavity 
depth  and  loading. 

A  second  moment  method  solution,  termed  a  hybrid  basis  solution, 
was  developed  specifically  for  wide  shallow  cavities  with  enough  loading 
to  support  at  most  one  surface  wave  mode.  The  hybrid  basis  solution 
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employs  entire  basis  functions  and  Galerkin's  method,  as  does  the  modal 
basis  solution,  but  expands  the  equivalent  magnetic  current  as  a  sum  of 
forced  wave,  two  surface  waves  (if  present)  and  several  transition 
waves.  The  transition  waves  are  exponential  functions  chosen  to  be  most 
significant  near  the  edges  of  the  aperture  in  the  surface  wave 
transition  region.  This  representation  was  found  to  yield  accurate 
results  for  the  hard  polarization  for  cavity  depths  <  .lOA  and  widths  > 
several  x  with  the  use  of  just  3-4  transition  waves.  Thus,  since  only 
about  10  unknowns  are  required  for  any  aperture  width  and  the  basis 
functions  are  not  significantly  more  complicated  than  are  the  modal 
basis  functions,  the  hybrid  basis  solution  is  an  extremely  efficient 
solution  for  wide  shallow  cavities.  For  the  soft  polarization  the 
surface  wave  mode  was  not  present  for  cavity  depths  <  .lOX  and  for  the 
choice  of  dielectric  material  considered.  Also,  the  transition  waves 
were  found  to  not  be  appropriate  for  this  polarization.  Thus,  only  the 
forced  wave  plus  two  exponentials,  chosen  so  that  the  forced  wave  would 
meet  the  proper  edge  condition,  made  up  the  hybrid  basis  functions. 
Fairly  accurate  results  were  obtained  for  cavity  depths  <  .10  X  and 
widths  >  several  x,  except  near  grazing  angles  of  incidence.  Again, 

hybrid  basis  solution  is  extremely  efficient. 

The  analytical  techniques  presented  here  should  prove  useful  in  the 
design  of  surface  wave  antennas  and  in  determining  the  RCS  contribution 
from  such  structures  as  well  as  from  loaded  and  unloaded  notch 
structures.  Future  research  could  include  modification  to  the  cavity 
Green’s  function  to  accomodate,  for  example,  impedance  side  walls. 
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Improvements  to  the  hybrid  basis  solution,  specifically,  may  be  made  by 
choosing  more  appropriate  transition  wave  basis  functions.  Such  basis 
functions  may  be  obtained  from  a  UTD  analysis  of  an  appropriate 
canonical  surface  wave  diffraction  problem,  or  perhaps  by  using  many 
pulses  or  piecewise  sinusoids  near  the  edges  of  the  aperture. 
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APPENDIX  A 


DERIVATION  OF  THE  INTEGRAL  EQUATIONS 


The  integral  Equations  (2.2)  and  (2.4)  for  the  region  I  (y>0)  and 
region  II  (y<0)  fields,  respectively,  may  be  derived  from  Green’s  second 
identity  for  scalar  fields  and  the  equivalence  theorem.  In  two 
dimensions  Green's  second  identity  has  the  form  [16] 


3  ^  3d; 

/  ('*'  In  -  yn)  =  //  ds  (A.l) 

L  s 


where  ij/  and  <(i  are  scalar  fields  and  s  is  a  surface  bounded  by  the 
contour  L.  Note  that  this  identity  is  a  statement  of  reciprocity  for 
the  two  fields. 


Let  4>  be  the  scalar  field,  Uj,  defined  by 

h 

E-  (P) 

Us  (p)  =  ^  ,  (A. 2) 

b  (p) 

and  let  be  the  appropriate  Green's  function,  G5(p,p').  To  derive 

h 

integral  equations  from  (A.l)  we  must  consider  the  wave  equations  and 
boundary  conditions  satisfied  by  u  and  G.  The  scalar  fields  satisfy 


jwp  J, 

(V2  +  k2)  Us  (p)  =  .  (A.3) 

h  juie 
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I 

where  J2,  are  impressed  sources  located  in  region  I  as  shown  in 
Figure  A.l.  On  perfect  conductors,  i.e.,  paths  L^.  and  C,  the  fields 
must  satisfy  the  boundary  conditions 

=  0  (A.4) 

L 

where  L  =  or  C  and  n  is  the  interior  normal  along  L.  The  Green's 
function  satisfies 

(V2  +  1(2)  Gs  (5.p')  =  -  6  (p.p')  ,  fA.5) 

h 

with  boundary  conditions  chosen  so  as  to  simplify  (A.l)  as  much  as 
possible. 

To  facilitate  determination  of  the  Green's  function  the  original 
problem  of  Figure  A.l  is  converted  into  an  equivalent  problem  by 
employing  the  equivalence  principle.  A  perfectly  conducting  sheet  is 
placed  in  the  aperture,  and  on  this  sheet  equivalent  source 
distributions  are  placed  on  either  side.  Then,  for  the  equivalent 
problem  in  region  I,  application  of  Green's  identity  along  with 

Equations  (A. 3),  (A.4)  and  (A. 5)  yields 

dUj  3G^ 

h  h 

/  (Gs  (p,p'}  —  (p')  -  Us  (?•)  7^  (p.p'))  dt' 

L  h  h 

d 

=  I!  Gs  (p.p')  1  '^0  }  ds'  +  Us  (p')  (A.6) 

Sj  h  joie  Mq  h 
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Figure  A.l.  Geometry  for  derivation  of  integral  equations. 


Choosing  the  boundary  conditions 


1  3Sh 

I  3n  =  0  on  Lg 


G5  becomes  the  half-space  Green's  function 
h 

Gs  (x,y;x'.y')  =  -j/4  (k  /U-x')2  +  (y-y')2) 

^  j/4  (k  y(x-x')2  +  (y+y')2) 


(A. 7) 


<A.8) 
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Assume  the  sources  in  region  I  to  be  z-directed  electric  and  magnetic 
line  sources  for  the  soft  and  hard  polarizations,  respectively,  located 
(xo»yo)*  Figure  A.l  and  Equations  (A. 7)  and  (A. 6), 

Equation  (A.fi)  becomes 


Uj  (p)  =  Us  (p)  +  Us  (p)  +  / 

h  h  h  -a 


3  _  _ 

Us  (p* )  3y '  fig  (p,p'  ) 
^^h  . 

(p')  (p.p') 


dx' 

y=0+ 


(A.q) 


where 


'Js  fp)  =  {  Hq  f  (k  *  y-Yo^^) 

iJs  (p)  =  (  Hq  }  (k^x-Xo)2  +  y+yQ)2) 


(A. 10a) 


(A. 10b) 


and  where  fis  given  by  Equation  (A. 8).  In  this  work  the  incident 
h 

field  is  assumed  to  he  (locally)  a  plane  wave,  thus  the  asymptotic  form 
of  Equation  (A. 10)  is  used.  This  completes  the  derivation  of  the 
integral  equation  for  the  fields  in  region  I. 

For  the  equivalent  problem  in  region  II,  application  of  Green's 
identity  and  Equations  (A. 4)  and  (A. 5)  (there  are  no  sources  in  this 
region)  yields 

3Us(p')  365 

_  _  h  h  _  _ 

/  ^Gs  (p,p')  “an  -  Us  (p')  ^  (p.p'))  dt' 
h  h 

O 

=  Us  (p)  .  (A. 11) 

h 
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Choosing  the  boundary  conditions  as  given  by  Equation  (A. 7)  becomes 
the  Green's  functions  for  the  enclosed  cavity  (derived  in  Appendix  B). 
Applying  the  boundary  conditions  Equation  (A. 11)  becomes 


a 

u^(p ' )  lyr-  6^  (p,p' ) 

>  I 

-a 

^  (p‘)  6h  (P.P’) 

(A.12) 


Where  G^  is  the  cavity  Green's  functions.  This  completes  the  derivation 
of  the  integral  equation  for  the  fields  in  region  II. 
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APPENDIX  B 


DERIVATION  OF  THE  CAVITY  GREEN'S  FUNCTIONS 


The  cavity  Green's  function  for  the  geometry  of  Figure  2.3  is 
easily  derived  by  the  method  of  separation  of  variables.  The  Green's 
functions  satisfy 

(|j^  Jp'  ^  ^r^  (x,y;x',y')  =  -6(x-x')  6(y-y')  (R.l) 

n 

with  the  boundary  conditions 


3Qh 

an 


=  0 


on  C 


(R.2) 


where  C  is  the  perimeter  of  the  cavity  defined  by  x=-a,  x=a,  y=-b,  y=0 

A 

and  n  is  the  interior  normal  on  C.  Consider  the  separation  of  variables 

[201 

Gs  (x,y;x',y')  =  K  Xj  (x,x')  •  Ys  (y,y')  ,  (R.3) 

h  h  h 

where  K  is  a  linear  operator  such  tha*  Equation  (R.3)  is  a  solution  to 
Equation  (B.l).  Substituting  (R.3)  into  (R.l)  suggests  the  ordinary 
differential  equations 


and  the  appropriate  boundary  conditions.  Consideration  of  the  form  of 
Equation  (R.7)  leads  to  the  choice  of  linear  operator 


(B.9a) 


(B.9b) 


(B.IO) 


Note  that  the  paths  of  integration,  and  Cy,  encircle  the 
singularities  of  X  and  Y,  respectively.  From  Equations  (B.7a)  and 
(B.9a),  Equation  (B.3)  can  be  shown  to  reduce  to 

G(x,y;x',y')  *  I  <^^(x)  ^^(x')  •  Y(y,y';X^)  .  (B.lla) 

m 


If,  instead,  Equations  (B.7b)  and  (B.9b)  are  used  the  result  is 

G(x,y;x',y')  I  (PnCy)  'I'pCy')  *  X(x,x’;X^)  .  (B.llb) 

n 

Consider  first  the  form  of  the  Green's  function  given  by  Equation 
(B.llb).  The  orthonormal  eigenfunctions  arising  from  Equation  (B.8b) 
with  the  boundary  conditions 

3'|'h,£  =0  on  y=-b,0  (B.12) 

_Ty  “ 

are 


s.t  (>> 

A 

tiry 

Sin  b 

h.t 

tny 

COS  b 
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where  l  is  an  integer  and  is  given  by  Equation  (2.13c).  In 
Equations  (B.llb)  X(x,x')  is  determined  from  Equations  (B.4a)  and  (B.6) 
with 

\  =  "b^  .  (R.14) 


The  solution  for  X(x,x')  is  constructed  in  two  intervals  and  the 
properties  of  the  Green's  function  at  x=x'  are  used  to  determine  the 

arbitrary  coefficients  involved.  That  is, 

IJ(x<;Xy)  •  T(x>;Xy) 

X(x.x';Xy)  = - ^  (B.15) 


where  U  and  T  are  independent  solutions  of 

(fe-  Mk2  I  ■  !>  (B-IM 

and  the  constant  W  is  (for  Equation  (B.16)) 

W(T,U)  =  T(x)  ir(x)  -  T'(x)  •  IJ(x)  .  (B.17) 


Solutions  to  Equation  (B.16)  are 


Us(’<)  1  _ 

sinT(a+x) 

1  ’  1 

sinY(a-x) 

(B.IR) 


for  the  soft  polarization  and 


Uh(x) 

coSY(a+x) 

Th(x) 

COSY (a-x) 

(B.ig) 


for  the  hard  polarization,  where  i  is  given  by  Equations  (2.13a)  and 
(2.13b).  The  constant  W  is 
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Ws  =  t  Y  sin2Ya  (B.20) 

h 

for  the  soft  and  hard  polarizations,  respectively.  From  Equations 
(B.18),  (B.19)  and  (B.20),  Equation  {B.15)  may  be  written  in  the  form 
-csc2Ya  , 

(x.x')  =  2y  [cosY(2a-|x-x' I )  +  cosy{x+x')]  .  (B.21) 


Finally,  (B.llb),  (B.13)  and  {B.21)  yield  the  results  given  by  Equations 
(2.12a)  and  (2.12b). 

Consider  now  the  form  of  the  Green's  function  specified  by 
Equation  (B.lla).  An  orthonormal  set  of  eigenfunctions  is  determined 
from  Equation  (B.8a)  with  appropriate  boundary  conditions.  Equation 
(B.4b)  is  then  solved  by  the  same  method  used  to  solve  Equation  (B.4a) 
in  the  previous  derivation.  The  final  results  for  the  Green's  functions 
are 

...  -  1  *  cscYb  .  . 

G5(x,y;x',y')  =  2a  I  sin  ^  (a-x)  sin  (a-x  ) 

4=1 

• [cosY(b-|y-y ' I )  -  cosY(b+y+y' )]  (B.22a) 


Gf,(x,y;x'  ,y' ) 


-  1  ®  CSCYb 

2a  I  1777  cos  2^  (a-x)  cos  JE 
4=0 


(a-x') 


•[cosY(b-|y-y' 1 )  +  cosY(b+y+y' )] 


(B.22b) 


where 


4  < 

4  > 


2ak 

TT 

2ak 

IT 


Re  (e^  u^) 

(B.23a) 

Re  (^p  Up)  . 

(B.23b) 
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APPEMDIX  C 


DERIVATION  OF 


In  this  appendix  we  derive  expressions  for  the  elements  of  the 
impedance  matrices  of  the  solutions  presented  in  this  thesis.  Consider 
first  the  hard  polarization.  From  Equations  (3.8)  and  (2.7),  let 

^  ^  H^^^(k|x-x' I )  -  Cp  (x.x')]  m(x’)  w(x)  dx'dx 


-a  -a 


(C.l) 


where  m=M^  and  w=Wj^  and  is  given  by  Equation  (2.12b).  Now  consider 
just  the  first  term  in  Equation  (C.l).  Let 


a  a 


^u  ■  /  /  (k I x-x' I )  m(x')  w(x)  dx'dx  .  (C.2) 

-a  -a 


This  expression  may  be  reduced  to  a  single  integral  over  the  Hankel 
function  through  the  coordinate  rotation  (the  so-called  Popovich 
transformation) 


x-x ' 

^  "  /2 


u’ 


X-t-X  ' 
/2 


(C.3a) 

(C.3b) 


The  resulting  integral  is 

/2a  /?a-lu|  (2)  u '  -u  u  '■♦•u 

=  /  /  j/2  (/2k|u|)  m[  w[  )  du'du  (C.4) 

-/zi  -/Za*\\i\ 
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From  Equations  (3.13)  and  (3.30)  we  see  that  the  basis  functions  and 


weighting  functions  are  simply  exponentials.  Thus,  let 
m(x')  =  eJ^^' 


(C.5a) 


w(x)  = 


(C.5b) 


Substituting  Equations  (C.5)  into  Equation  (C.4)  and  evaluating  the 


integral  over  u', 

a  /px  sin[(f+g)(a-u)] 

^u  ^  ^(2ku)  cos(f-g)u - - du  . 


(C.6) 


For  ku  near  zero  the  small  argument  form  of  the  Hankel  function  [17], 


Hq^Vp)  -  1  -  j  r  [^)  .  kp  >  0 


(C.7) 


where  y  =  1.781,  is  used  so  that  the  singularity  at  the  origin  may  be 
integrated  out.  Substituting  (C.7)  into  (C.6)  for  u  <  .05  yields  the 
approximate  result, 

Au. 

Au  ^  sin[(f+g)(a--2)  ^  2  ,YkAu,, 

^u  Q  *  4jcos[(f-g)  2]  •  •  [1-j  7  tn  [“ i")]  • 

(C.7.  ) 

For  ku  large  (2ku>5)  the  large  argument  form  of  the  Hankel  function,  as 
given  by  Equation  (2.8.1),  is  used  to  express  most  of  Equation  (C.6)  in 
terms  of  the  complex  error  function  [17], 


2  -t2 


(C.8) 
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by  employing  the  relations 


.-YX 


/  dx  =  /tTTy*  erf  Ax' 

_ erf  Ax  /x* 


-tx 


/  AT  e'^^dx  =  /n/Y 


2y 


(C.9a) 

(C.9b) 


For  f+g  +  0  the  relevant  portion  of  Equation  (C.6)  becomes 


0.4 


f+g  v  irk 


J(i^+g)a 


/*  ^-7=-  •  *  e  -^29")  du 

C.4  ^ 


■j2ku 


j‘  .  (^j2fu  ,  jJ29u, 

0.4 


-j2ku 


(C.IO) 


which  is  easily  evaluated  with  Equation  (C.9a).  For  f+g=0  the 
corresponding  result  is 


g-j2ku 
(a-u)  -7^ 


(  ^  e-j(f-g)u)  ^ 

(C.ll) 


which  is  easily  evaluated  with  Equations  (C.9a)  and  (C.9b).  The 
remaining  portion  of  the  integration  in  Equation  (C.6),  from  u=.05  to 
u=0.40,  is  evaluated  numerically  by  the  trapezoidal  rule. 

Efficient  code  for  computation  of  the  Hankel  function  and  complex 
error  function,  as  well  as  an  efficient  program  for  complex  matrix 
inversion,  were  provided  by  Prof.  J.H.  Richmond.  These  routines  are 
discussed  in  [12]. 
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Next  we  consider  the  second  term  in  Equation  (C.I.).  denoted  Zl, 


where 


-e  /  /  G.  (x,x')  m(x')  w(x)  dx'dx 

-a  -a 


(C.12) 


Substituting  Equations  (C.5a),  (C.5h)  and  {?.12b)  evaluated  at  y=y'=n 
into  Equation  (C.12),  the  resulting  double  integral  may  be  evaluated  in 
a  straight  forward  manner.  The  resulting  sum  is 
Ef.  sin(f+g)a  -  I  1 


4=0  Y^-f^ 


•  csc2Ya  1  1“  sin(Y+g)a 

■  *  ‘r  ,L  ^ 


sin(Y-g)a-| 
sin(Y-f)a  _! 


{C.13) 


where  y  is  a  function  of  l  and  is  given  by  Equation  (2.13).  The  first 
sum  in  this  expression  may  he  put  in  closed  form  through  the  use  of  the 
relation  [22] 


«  1  1  It 

Jl  *  2^  ■  il 


(C.14) 


Thus,  the  first  sum  in  (C.13)  becomes 

»  1  1  b  cot{b  /  k2epUp-f2) 

^lo  *  Y^-f^  ~  ^  /'k^e^Up-f^ 


(C.15) 


For  the  cavity  depths  of  interest  in  this  work  {<  1/4  wavelength)  only  a 
few  terms  of  the  second  sum  in  Equation  (C,13)  are  needed  for 
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convergence,  regardless  of  cavity  width.  Thus,  with  the  use  of  Equation 
(C.15)  a  highly  convergent  formulation  of  Zl  has  been  achieved. 

An  equivalent  expression  for  may  be  obtained  by  using  the 
alternative  form  of  the  cavity  Green's  function  given  by  Equation 
(B.22b).  Substituting  this  expression  evaluated  at  y=y'=0,  along  with 
Equations  (C.5a)  and  (C.5b),  into  (C.12)  yields  the  result 


cotYb 


sin(f+p)a 
f+P  ' 


sin(g-t-p)a 

g+p 


sin(g-p)a  sin(f-p)a  sin(g-t-p)a  £  sin(g-p)a 

+  g-p  _  f-p  _  g+p  g-p 


(C.16) 


where  y  is  given  by  Equation  (B.23)  and 
An 

P  =  ?a 


(C.17) 


No  convenient  method  was  found  for  reducing  any  portion  of  the  sum  in 
Equation  (C.16)  to  closed  form,  making  this  representation  of  Zl  far 
less  convergent.  The  number  of  terms  required  in  the  summation  is 
directly  proportional  to  the  width,  and  is  typically  equal  to  10-15 
times  the  width  in  wavelengths. 

Now  consider  the  soft  polarization.  In  view  of  Equations  (3.8)  and 


(2.8)  let 


i  f‘ 

=  2  J 


-a 


32 

TP 


Hi^V/(x-x')2+y"27 


m(x')  w(x)  dx'dx 
y  =  O'*- 


(C.18) 
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where  m(x*)  and  w(x)  are  basis  and  weighting  functions  such  that 


m  (±a) 
w  (ta) 


=  0 


(C.19a) 

(C.19b) 


From  the  differential  equation  which  the  2-D  half-space  Green's  function 
satisfies,  and  noting  that  y=^y',  we  obtain  the  relation 

Ip-  (kl^p'l)  -  -(Ip  +  k2)  (kip.p'l)  (C.20) 

where  |p-p'|  Is  given  by  Equation  (2.3).  Substituting  this  into 
Equation  (C.18)  we  then  integrate  by  parts  to  remove  the  derivatives  on 
the  Hankel  function.  The  end-point  terms  reslting  from  this  integration 
can  be  shown  to  vanish  from  a  consideration  of  the  ege  condition  and  the 
small  argument  of  the  Hankel  function,  as  given  by  Equations  (2.14)  and 
(C.7),  respectively,  and  from  the  relation 

lim  x^  tnx  =  0  ,  t  >  0  .  (C.21) 

x*0 

The  resulting  form  of  Z^  Is 


~  ^  f  f  H^^^(k |x-x' I )  m(x')  w(xl  dx'dx 

-a  -a 


+  ■?  /  /  H^^'{k|x-x' !)  TP  (x* )  (x)  dx'dx  .  (C.22) 


-a  -a 


From  Equations  (3.15)  and  (3.22)  It  Is  apparent  that  the  basis  and 
weighting  functions  are  a  sum  of  exponentials.  Thus,  the  evaluation  of 
Zy  for  the  soft  polarization  Is  fundamentally  the  same  as  for  the  hard 
polarization. 
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Finally,  consider  the  second  term  in  Equation  (2.8).  Let 


1  a  a  32 

=  ur  I  /  [ay^  G^(x,y;x'.y')]  m(x' )w(x)dx'dx 

y'=0- 
y  =0" 


(C.23) 


where  m(x')  and  w(x)  satisfy  Equation  (C.19).  Substituting  Equation 


(2.2)  into  (C.23)  yields 


it2  - 

^  ih 


csc2ya 


^  tiry* 

42  cos  "h  cos  •  I 


(C.24) 


where 


I  =  /  /  [coSY(2a-lx-x' I )  -  cosy(x+x*)]  m(x')w(x)  dx'dx  (C.25) 

-a  -a 

and  Y  is  a  function  of  l  as  given  by  Equation  (2.13).  To  obtain  a 
convergent  summation  Equation  (C.25)  must  be  integrated  by  parts,  which 
yields  the  result 
1  39 

^  ~  I  I  [coSY(2a-|x-x' I )  +  cosy(x+x' )]  “35^  3^^  dx’dx 

-a  -a 


2sin2Ya  ® 


/  m(x)w(x)  dx 

-a 


(C.26) 


Reecalling  that  m(x')  and  w(x)  are  a  sum  of  exponentials  for  this 
polarization,  we  substitute  Equations  (C.5)  into  Equation  (C.26)  and 
evaluate  the  integrals  in  a  straightforward  manner.  Employing  the 
rel ation 


r  1  4iiy  tiry' 

1  IZTb  cos  h  cos  h  =  «(y-y') 
1=0 


(C.27) 
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along  with  Equation  (C.14),  and  noting  that  y+y'.  Equation  (C.24) 
becomes 


csc2Ya 


“  1 


r  •  /  sin(Y^^g)a  .  sin(Y-g)a 

•  [sin(Y+f)a  Y+g  “  sin(Y-f)a  y.g  ] 


(C.28) 


where 


Z,  = 


2  sin(f-*-g)a 

Ur  f+g 


A^e  u 
r  r 


•  cot  b  /k^e^Up-f^  -  (1  +  f"  A^EpU^  cot  b  A^EpUp) 


(C.29) 


A  second  form  for  Zl  is  obtained  by  employing  the  alternative  form 
of  the  cavity  Green's  function  given  by  Equation  (B.22b).  Substituting 
this  expression  into  Equation  (C.23)  and  noting  that  y4=y',  Zl  becomes 
-1  * 

Z.  =  YCOtYb  •  I  ,  (C.30) 

t=l 


where 

a  a 

I  =  /  /  sin  p(a-x')  sin  p{a-x)  m(x')  w(x)  dx'dx  (C.31) 

-a  -a 


and  p  is  given  by  Equation  (C.17}.  Integrating  by  parts,  {C.31)  becomes 

1  a  3  ,  ,  3m(x')  3w(x)  . 

I  =  ^  /  /  cos  p(a-x')  cos  p(a-x)  -gT*  ^  dx  dx  .  (C.32) 

-a  -a 

Substituting  Equations  (C.5)  into  (C.33)  and  evaluating  the  double 
integrals,  Equation  (C.32)  is  substituted  into  (C.30)  and  Zl  becomes 
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fg  ?  YCOtYb 

^  Ur  a 


sin(f+p)a 

f+p 


sin(g■^p)a 
(-1)  g+p 


sin(g-p)a~ 

g-p 


sin(f-p)a 

^~f-P 


'sin(g-t-p)a  .  s1n(g~p)a 

g+p  \“1)  a-D 


g-p 


(C.33) 


where  y  is  given  by  Equation  (B.23)  and  p  is  given  by  Equation  (C.17). 

Equations  (C.28)  and  (C.33)  possess  the  same  convergance  properties 
that  the  corresponding  hard  polarization  equations  have.  Thus, 

Equations  (C.28)  and  (C.29)  yield  a  far  more  convergent  representation 
than  does  Equation  (C.33)  for  wide  cavities. 
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